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ABSTRACT 

We present an analysis of the evolution of the central mass-density profile of massive elliptical 
galaxies from the SLACS and BELLS strong gravitational lens samples over the redshift interval 
z w 0.1-0.6, based on the combination of strong-lensing aperture mass and stellar velocity-dispersion 
constraints. We find a significant trend towards steeper mass profiles (parameterized by the power-law 
density model with p oc r~ 7 ) at later cosmic times, with magnitude d (7) jdz — — 0.60±0.15. We show 
that the combined lens-galaxy sample is consistent with a non-evolving distribution of stellar velocity 
dispersions. Considering possible additional dependence of (7) on lens-galaxy stellar mass, effective 
radius, and Sersic index, we find marginal evidence for shallower mass profiles at higher masses and 
larger sizes, but with a significance that is sub-dominant to the redshift dependence. Using the results 
of published Monte Carlo simulations of spectroscopic lens surveys, we verify that our mass-profile 
evolution result cannot be explained by lensing selection biases as a function of redshift. Interpreted 
as a true evolutionary signal, our result suggests that major dry mergers involving off-axis trajectories 
play a significant role in the evolution of the average mass-density structure of massive early-type 
galaxies over the past 6 Gyr. We also consider an alternative non-evolutionary hypothesis based on 
variations in the strong-lensing measurement aperture with redshift, which would imply the detection 
of an "inflection zone" marking the transition between the baryon-dominated and dark-matter halo- 
dominated regions of the lens galaxies. Further observations of the combined SLACS+BELLS sample 
can constrain this picture more precisely, and enable a more detailed investigation of the multivariate 
dependences of galaxy mass structure across cosmic time. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: structure 
— gravitational lensing: strong 



1. INTRODUCTION 

Massive elliptical galaxies play a starring role in both 
galaxy evolution and cosmology. In the former con- 
text, they are the end-products of hierarchical galaxy 
merging Ce.g.. iToomre fc Toomrelll972t ISchweizerlllQSa 



IWhite fc F rcnk 1991; Kauff mann et~aT1 [l993: Co le et all 
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I2000T ) . and in the latter context they are highly biased 
tracers of large-scale structure that are easily visible 
at cosmological dis tances (e.g., lEisenstein et al.l 120051 : 
iPercival et al.ll2007|) . Hence cosmological spectroscopic 
survey samples such as the lu minous red galaxy sample 
(LRG: lEisenstein et al.1 12001ft of the Sloan Digital Sky 
Survey ('SPSS: lYork et al.l 12000ft and the Baryon Os- 
cillation Spectroscopic Survey of the SPS S-Ill (BOSS: 
lEisenstein et al.l 120111 : iPawson et al.1 12012ft provide un- 
cqualcd opportunities to study massive galaxy evolution. 

Important problems in massive galaxy evolution re- 
main open, including quant ifying the role of merg- 
ers in their evolution (e.g.. Ivan Pokkum et al.1 119991 : 
iKhochfar fc Burkerti2003t iBell et al.ll2006| ). deducing the 
histor y of their ste l lar p o pulations over cosmic time 
(e.g., (Thomas et all 120051 ; iMaraston et~aT] 120091) . and 
identifying the m echanism for thei r observed evolu- 
tion in size (e.g lPaddi et al.l 120051: iZirm etal l 120071 : 



Ivan der Wel"et^[200l " lvan Pokkum et al.l 120081 ). The 

persistence of these and other problems is due in large 
part to the difficulty of making precise galaxy mass 
measurements at cosmological distances. Stellar mass 
measurements based upon integrated photometric and 
spectroscopic diagnostics, in addition to being insen- 
sitive to dark- matter content, are limited by uncer- 
tainties in the initial mass function and othe r stellar- 
population parameters (e.g.. lConrov et ai1l2009l) . and de- 
tailed dynamical modeling methods applicable to early- 
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type galaxies in the lo cal universe (e.g.. iGerhard et all 
[200lHCaDPellari et al.ll2006h suffer from mass-profile and 
orbital-anisotropy degeneracies at cosmological distances 
where higher-order moments of the line-of-sight velocity 
profile can no longer be measured reliably. 

The phenomenon of strong gravitational lensing per- 
mits mass measurements in the central regions of early- 
type galaxies to the few-percent level. Either on 
its own (e.g . iWarren fc Dvel 12001 iRusin et al.l 12001 
IWavth et all 12005ft or in combination with stellar- 
dynamical constraints (e.g.. jTreu fc Koopmansl 120021 
120041 : IKoopmans k, Treul 120031) . strong lensing provides 
the most powerful tool for the measurement of galaxy 
masses and density distributions across cosmic time. 
The use of strong lensing for the study of galaxy struc- 
ture and evolution has historically been constrained by 
the availability of significant lens samples. By iden- 
tifying gravitational lens candidates spectroscopically 
from wit hin the SPSS database, t h e Sloan Len s ACS 
(SLACS: IBolton et all [20061 120081: lAuger et al.l 120091) 
Survey has identified the largest single sample of con- 
firmed strong-lens galaxies, but the small redshift range 
within the SLACS sample alone limits its utility for the 
study of galaxy evolution. An initial combination of the 
SLACS sample wi t h lens es at significantly higher redshift 
(|Koopmans et all |2006|) detected no significant mass- 
structure evolution, but used a relatively small and het- 
erogeneous higher- redshift sample. Recently, the CFHT 
Strong Lens Legacy Survey (SL2S: iCabanac et all 120071 : 
iMore et all 120111 : iGavazzi et al.ll2012[ ) provided a larger 
sample of strong lenses at higher redshift, which yielded 
a tentative detection of structur al evolution whe n com- 
bined with the SLACS sample (jRuff et all I20ITI) . The 
SL2S sample is most fundamentally limited by a lack of 
source-galaxy redshifts for many systems, compromising 
the precision with which angular lensing observables can 
be translated into physical mass constraints. 

In this paper, we present a galaxy mass-density profile 
evolution analysis that combines the SLACS lens sample 
with lenses recently discov ered by the BOSS Emi ssion- 
Line Lens Survey (BELLS: iBrownstein et al1l2012L here- 
after Paper I). The BELLS lenses are selected from 
the BOSS survey with the same spectroscopic selection 
method employed by SLACS, and are comparable in stel- 
lar mass to the SLACS lens galaxies (see £|4.2l below.) For 
both SLACS and BELLS samples, spectroscopic lens and 
source galaxy redshifts are available for all systems, and 
SDSS/BOSS spectra provide stellar velocity-dispersion 
data that may be combined with lensing measurements 
to constrain the lens-galaxy mass-density profiles. 

This paper is organized as follows. Section [2] gives an 
overview of the observational data that we employ. Sec- 
tion [3] describes our inference methodology and presents 
our basic evolutionary measurement result. Section [4] 
presents an investigation of possible systematics and de- 
pendences beyond redshift evolution, includ ing those as- 
sociated with stellar velocity dispersions ( 34.1[) . photo- 
metric parameters ( £!4.2j) . strong- lensing measurement 
apertures ( §4.3[) . and lensing selection effects ( §4.4[) . Fi- 
nally, $5] provides a discussion and conclusions. We as- 
sume throughout a general-relativistic FRW cosmology 
with parameters Hq = 70kms _1 Mpc -1 , fiji = 0.3, and 
SI A = 0.7, and note that our results are only very weakly 
sensitive to the exact values of these parameters. 



2. OBSERVATIONAL SAMPLES 

The SLACS and BELLS lenses were selected from 
imaging and spe ctroscopic data co l lected at the 2.5-m 
SDSS telescope (|Gunn et all fl998l 12006ft . all of which 
are included in the public SDSS-HI Data Release 9 
(|Ahn et al.l 120121) . For the high-resolution imaging 
follow-up necessary for strong-lcnsing analysis, all sys- 
tems in this work were also observed through the F814W 
(/-band) filter using the Wide-Field Channel (WFC) of 
the Advanced Camera for Surveys (ACS) aboard the 
Hubble Space Telescope (HST) . For our SLACS data 
set, we use the 57 early-type grade-A lenses with lens- 
ing mass models presented in IBolton et al.l (|2008f ). We 
do not in clude data fo r the a dditional SLACS lenses re- 
ported in lAuger et al.l (|2009( ) , since that subset was ob- 
served only with the WFPC2, and we wish to maintain 
as much uniformity as possible. Since the uncertainty in 
our evolutionary analysis is primarily set by the BELLS 
sample at the higher redshift end, the omission of these 
additional SLACS systems has little statistical effect. For 
our BELLS data set, we use the 22 early-type grade-A 
lenses presented in Paper I. For all lenses, we use SDSS 
(for SLACS) and BOSS (for BELLS) spectroscopic data 
to derive stellar velocity-dispersion likelihood functions. 
For 6 of our BELLS targets, multiple independent spec- 
troscopic observations are available. We combine the 
velocity-dispersion likelihood information from these re- 
peat observations in our analysis below. 

3. BASIC EVOLUTION MEASUREMENT 



Following IKoopmans etldl (|2006l I2009D . we constrain 
the central mass-density profile of the lens galaxies via 
a Jeans-equation analysis constrained by both strong- 
lensing aperture masses and stellar kinematics. In this 
analysis, we model all lenses using a stellar luminosity 
profile embedded in a total mass-density profile (i.e., 
stars plus dark matter) parameterized as p cx r -7 . 

For the parameterization of the lu minosity profile, w e 
use the "Nukcr" profile described by lLauer et all (|1995l ) , 
which is a broken power-law with a transition of vari- 
able softness between inner and outer regions. We fit the 
PSF-convolved Nuker profile directly to the HST J-band 
images inside a circular region of radius 6", and optimize 
the parameters non-linearly until convergence. This cen- 
tral region is chosen so as to focus on accurately matching 
the luminosity-profile model to the imaging data over the 
approximate range of radii probed directly by our lens- 
ing and kinematic observables. This approach allows a 
more accurate Jeans-equa tion analysis than one based 
on the globally optimized Ide Vaucouleursl ( 11948 1) mod- 



els which were used f or the measur e ment of magnitudes 
and effective radii in IBolton et al.l (|2008h and Paper I. 
However, we verify further below that our results do not 
depend significantly upon our choice of parameterized 
l uminosity-profil e form . 

I Graham et al.l (120031) demonstrate that the parameters 
of the ISersid (jl968h ~nrofile are more robust and physi- 
cally meaningful than those of the Nuker profile for the 
global characterization of galaxy surface-brightness pro- 
files. We emphasize here that the role of the Nuker model 
in our mass-profile analysis is simply to provide a flexi- 
ble, convenient, and sufficiently accurate parameterized 
mathematical model for the luminous tracer distribution 
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Fig. 1. — Comparison of \ 2 goodness-of-fit statistics for Sersic 
and Nuker parameterized models to SLACS and BELLS lens galax- 
ies used in this work, as fitted within circular regions of radius 6" 
centered on the lens galaxies. 

over the local range of radii probed by our strong-lcnsing 
and stellar-dynamics observations. We make no physi- 
cal interpretation of the fitted Nuker model parameters, 
which were originally conceived to describe the core re- 
gions of galaxies. Figure Q] shows a comparison of the 
X 2 goodness-of-fit statistic for Sersic and Nuker fits to 
the inner regions of our SLACS and BELLS lens galax- 
ies. The quality of the Nuker fit is as good as or better 
than that of the Sersic fit in every case, with the typical 
improvements in x 2 being highly significant Q Hence we 
adopt the Nuker profile as our preferred parameterization 
within the lensing-dynamical analysis, but we return to 
the Sersic model to characterize the distribution of lens- 
galaxy global profile shapes in £14.21 below. 

For each galaxy, we take the measured Einstein ra- 
dius as a projected aperture-mass constraint on the to- 
tal mass-density profile. We deproject this mass con- 
straint using the power-law mass model for a gridded 
range of profile values 1.1 < 7 < 2.9. Assuming a con- 
stant and isotropic velocity-dispersion tensor, we apply 
the deprojected mass constraint and solve the spheri- 
cal Jeans equation at each point in the grid of 7 val- 
ues for the radial (squared) velocity-dispersion profile of 
a tracer population given by the Abel-deprojected lu- 
minosity profile for the lens galaxy. We then reproject 
the luminosity-weighted linc-of-sight component of this 
profile back into two dimensions to predict the line-of- 
sight velocity dispersion profile across the lens galaxy. 
The luminosity-weighted squarcd-dispcrsion is then in- 
tegrated over an aperture corresponding to the SDSS 
(3" diameter) or BOSS (2" diameter) spectroscopic fiber 
aperture blurred by 1'.'8 FWHM seeing to predict (74(7), 
the velocity dispersion that would be measured by SDSS 
or BOSS for lens i given a logarithmic mass-density slope 
parameter 7. 



We revisit the method of measurement of the stel- 
lar velocity dispersion a from SDSS and BOSS spectra, 
in order to (1) ensure the greatest uniformity between 
the treatment of the SLACS and BELLS samples, (2) 
maximize the signal-to- noise ratio (S/N) of the BOSS 
velocity-dispersion measurements by using the full range 
of spectral data, and (3) propagate the full velocity- 
dispersion likelihood information in our analysis. This 
is well motivated since velocity-dispersion uncertainties 
dominate the statistical error budget, particularly for 
the lower signal-to-noise ratio (S/N) BOSS spectra. To 
do this, we develop a new velocity-dispersion measure- 
ment code that explicitly projects high-resolution stellar 
templates through the wavelength-dependent line-spread 
function recorded on a fibcr-by-fibcr basis for all BOSS 
and SDSS spectra. We also derive a new set of stella r 
templates from the Indo-US library (jValdes et al.ll200"4T ) . 
which we furthermore patch (over telluric bands) and 
extend (to redder and bluer wavelengths) using best-fit 
model atmospheres se lected from the POLLUX database 
(|Palacios et al.l l2010f ). This stellar tem plate genera- 
tion p rocedure is described in detail in iBolton et al.l 
(|2012( ). and is used primarily to generate stellar clas- 
sification templates for use in the BOSS spectroscopic 
pipeline. For the velocity-dispersion analysis, we only 
use stars of spectral type A through K. We perform 
a principal-component analysis decomposition of these 
templates, and retain the top 5 eigenspectra as our 
velocity-dispersion template basis. At each trial ve- 
locity dispersion, we fit the lens-galaxy spectrum as a 
linear combination of broadened basis templates (thus 
marginalizing to some extent over stellar-population un- 
certainties) plus a quartic polynomial, and accumulate 
the function xf(<7i) for each lens galaxy i. We make 
use of these full xf functions, rather than point esti- 
mates of velo city disper s ion, in our mass-profile analysis. 
As shown in iShu et aTl ()2012[ ) , this approach allows us 
to combine data from multiple galaxies and recover un- 
biased estimates of population distribution parameters 
even in the limit of low S/N in individual spectraFl 

To derive our mass-profile constraints, we compose the 
dynamical information encoded by Xi( a i) with the lens- 
ing information encoded by 0^(7) to obtain the proba- 
bilities for the spectroscopic data of individual lenses 
given a value for the logarithmic mass-profile slope 7: 



Pi(di| 7 ) oc exp[-Xi(<7i(7))/2]. 



(1) 



To model the distribution of 7 values within the popu- 
lation, we parameterize the conditional probability den- 
sity function (pdf) of 7 at a given lens rcdshift by a Gaus- 
sian whose mean value evolves linearly with redshift: 



P{l\z; 70, lz, S-y) 



x exp 



[7-(7o + 7s(*-0-25))] s 
24 

Here, 70 is the mean 7 value at redshift ; 



(2) 



0.25, 



11 Sec Byun c t a l, 1996 for a formula expressing the Sersic profile 
as a limit of the Nuker profile. 



12 Although we do not use point estimates of velocity dispersion 
in our analysis, we can quote a median fractional error in velocity 
dispersion of about 12% and a median absolute error of about 
22kms — 1 for the BELLS lens sample, defined by the half- width of 
the A\ 2 = 1 interval about the minimum-x 2 value. 
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is the evolution of the mean 7 with rcdshift, and s 7 is 
the intrinsic scatter in 7 at all redshifts. The choice 
to parameterize the redshift dependence relative to the 
intermediate value of z = 0.25 is made so as to minimize 
covariance between the model parameters 70 and j z . 

With these ingredients, we can now express the prob- 
ability of the observed spectroscopic data for lens galaxy 
i given the evolution-relation parameters as 

p(di\7o,7z,s 7 ) = J d-fp l {d l \-f)p(-f\z i ;-f ,j z ,s 1 ), (3) 

with the two factors in the integral given by Equations [1] 
and [2] Finally, the likelihood function for the population 
parameters is obtained from the product over all lens 
galaxies in the sample: 

£{lo,lz, s 7 \{di}) =p({di}|7o,7 Z! s 7 ) 

= JJp(d,|7o,7 2 ,s 7 ), (4) 

i 

with the individual factors obtained from Equation [3] 
We grid the space of (70, 7z, s 7 ) and calculate Equation^] 
at each grid point, thereby mapping out the likelihood 
surface for our evolution-relation parameters. We assume 
a prior p(7o,7 z ,s 7 ) that is uniform in 70, j z , and logs 7 
to define the posterior probability 

p(7o,7z,S7l{ d *}) K £ho,lz,s~ f \{d l })p(-y ,^ z ,s~ f ). ( 5 ) 

Figure [2] shows maximum-likelihood estimates of 7 
for the individual lenses (for illustrative purposes), as 
well as the posterior probability contours for the zero- 
point 70 and evolution 7 Z of the population mean mass- 
density slope for both the Nuker profile-based calcula- 
tion and the alternative de Vaucouleurs profile-based cal- 
culation. Characterizing the one-dimensional measure- 
ments via marginalization and Gaussian fitting, we find 
7 2 = -0.60 ± 0.15, 70 = 2.11 ± 0.02, and log 10 (s 7 ) = 
log 10 (0.14) ± 0.06dcx. Thus, we find a slightly super- 
isothermal avera ge pro file (as in iKoopmans et al.l 120091 
and lAuger et al.ll2010f l. and a very significant signal of 
evolution in the total mass-density profile of strong-lens 
galaxies in the sense of a "steeper" (more centrally con- 
centrated) profile with increasing cosmic time (decreas- 
ing redshift). 

To verify the robustness of our result against changes 
in the details of our analysis procedure, we run the 
same mass-profile evolution calculation using the best-fit 
de Vaucouleurs models to describe the luminosity profiles 
of the lens galaxies, and also using Nuker models fitted to 
the HST imaging data with flat rather than noise-based 
pixel weighting. The resulting parameters and associated 
uncertainties for these alternate calculations are given in 
Table [U and show only insignificant differences with re- 
spect to our reference model results. We also verify the 
degree of information added by the BELLS sample by 
repeating the mass-density profile evolution analysis us- 
ing the SLACS sample alone. This computation yields 
a value for the evolution of the mean logarithmic mass- 
density profile slope of 7 Z = —0.61 ± 0.26: consistent 
with the SLACS+BELLS result, but at much lower sig- 
nificance given the reduced rcdshift baseline. The other 
parameters of this SLACS-only analysis are also listed in 
Table [1] 



4. INVESTIGATING NON-EVOLUTIONARY EFFECTS 

Before discussing our basic result further, we explore 
whether the mass-density profile redshift dependence 
that we measure could reflect the signature of some other 
physical dependence that is introduced through sys- 
tematic differences between the higher-redshift (BELLS) 
and lower-redshift (SLACS) components of our sample. 
For reference, Figure [3] shows the variation of multi- 
ple other parameters with redshift within the combined 
SLACS+BELLS lens sample that is used in this work. 

4.1. Velocity- dispersion effects 

The first alternative that we examine is a systematic 
variation of the typical sample velocity dispersion with 
redshift. Since the statistical velocity-dispersion uncer- 
tainties and the mass-profile uncertainties are almost per- 
fectly correlated, we do not attempt to fit and interpret 
7 versus cr, but rather we examine the variation of a with 
redshift. First, we correct the X 2 ( a ) velocity-dispersion 
baselines from the observational aperture velocity disper- 
sion a for each lens i to the estimated value er within 
an apertu re of one effective radiu s, using the empirical 
formula of iCappellari et ahl ()2006f ): 

o- s ,i = (R ohs /R c ) om a t , (6) 

where R obs = l'/5 for SDSS and l'.'O for BOSS. The ap- 
plied correction (i? b s /-R e ) ' 06 — 1 within the sample has 
a (signed) median of —1% and an RMS absolute value 
of 3%, since the effective radii of the lenses are gener- 
ally comparable to the angular radius subtended by the 
spectroscopic fiber. Although the correction relation is 
derived from galaxies in the local universe, the correc- 
tions are small enough that we assume any higher-order 
redshift-dependent effects to be most likely negligible. 

We fit for the distribution of a e values within the sam- 
ple as a function of redshift using the same method as for 
the 7 distribution above. We consider an evolving log- 
normal model for the population distribution in analogy 
to Equation [2] as follows: 

P( lo Sio °"c|z; too, m z ,Sa) = 

\/2irScr 

f [log 10 <7„- {m + m z (z-0.25))} 2 \ 
xex p \ 2 ^ j, {() 

for er c in kms -1 . Our notation follows IShu et all (|2012[ ) 
in using m to denote the peak in the population pdf 
of log 10 tr e values, with mo representing the value at 
z = 0.25 and m z parameterizing its evolution with red- 
shift. We constrain the parameters of this model as de- 
scribed via Equations [TH5] above, using xf (c e ) instead of 
X?(ci(7)). We find best-fit values of m = 2.39 ± 0.01, 
m z = -0.033 ± 0.053, and log 10 s a = log 10 (0.06) ± 
0.04 dex. Hence, we find no significant evidence of evo- 
lution in the mean value of the distribution of a c values 
within our sample of lenses. 

The upper left panel of Figure [3] shows the maximum- 
likelihood cr estimates for the SLACS and BELLS early- 
type lens galaxies as a function of redshift, along with the 
best-fit evolution relation and intrinsic scatter. While 
there is no significant evolution overall, we do see several 
possible outliers. To gauge the effect of these outliers 
on our 7 evolution analysis, we re-do our mass-profile 
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Fig. 2. — Left: Minimum-^ 2 values for the logarithmic total mass-density profile slope 7 for SLACS (black diamonds) and BELLS (blue 
squares) lenses. Error bars indicate A\ 2 = 1. The solid line shows the best-fit relation, gray lines indicate the "1-sigma" error in the slope 
and zero-point of this relation, and dashed lines indicate the best-fit intrinsic RMS population scatter. Red crosses indicate systems with 
a maximum-likelihood log-velocity dispersion (log 10 a e , in kms~ 1 ) cither greater than 2.5 or less than 2.2 (see Figurc[3]) Data points and 
error bars are for illustrative purposes only: the population parameter fits are done using the full x 2 (l) function for each lens, as described 
in ij3] Right: Posterior probability contours enclosing credible regions for the zero-point and evolution of the logarithmic mass-density 
slope parameter 7. Black curves are for the Nuker profile-based analysis; gray curves are for a de Vaucouleurs profile-based analysis. 



TABLE 1 

Mass-profile evolution parameters under various fitting scenarios 



Subset 70 7z log 10 cr 7 



All, statistical pixel weight 


2.11 


± 


0.02 


-0.60 


±0.15 


lot 


III 


(0.14) ±0.06 


dex 


All, flat pixel weight 


2.12 


± 


0.02 


-0.61 


±0.15 


lot 


ill) 


(0.14) ±0.05 


dex 


De Vaucouleurs image model 


2.13 


± 


0.02 


-0.58 


±0.15 


lot 


ill) 


(0.15) ±0.05 


dex 


SLACS lenses only 


2.11 


± 


0.02 


-0.61 


±0.26 


lot 


ill) 


(0.13) ±0.06 


dex 


Velocity-dispersion errors < 15% 


2.12 


± 


0.02 


-0.58 


±0.15 


lot 


S10 


(0.14) ±0.06 


dex 


2.2 < log ln fffkms" 1 ) < 2.5 


2.11 


± 


0.02 


-0.58 


±0.15 


l°f 


Sin 


(0.14) ±0.06 


dex 



evolution computation excluding the four lenses whose 
maximum-likelihood log 10 a e estimates (in kms -1 ) are 
either less than 2.2 or greater than 2.5 (corresponding to 
158 kms -1 and 316 kms -1 respectively). In this case, we 
find a value for the evolution of the population mean 7 
of 7 2 = —0.58 ± 0.15, which differs only insignificantly 
from the value we find when including all lenses. Table [T] 
contains the full parameter set from this fit. 

Changes in the redshift and spectroscopic S/N of the 
spectra from which the velocity-dispersions are measured 
are not likely to contribute to our observed signal. If 
such effects were driving our evolutionary measurement, 
we would expect to see them reflected in a significant ap- 
parent evolution in the velocity-dispersion distribution of 
our combined sample with redshift. However, as shown 
above, we do not see this effect to any significant degree. 
Furthermore, the hierarchical likelihood-stacking method 
that we use for incorporating velocity-dispersion infor- 
mation into our cur rent popu l ation analysis has been 
tested and verified in lShu et al.l (|2012fl . Specifically, that 
work has shown that increasing the redshift and simulta- 
neously decreasing the S/N of the spectra used for kine- 
matic analysis does not introduce any bias into the esti- 
mated parameters of the population velocity-dispersion 
function estimated from multiple spectra. If we exclude 
the 4 out of 79 lenses within our sample whose esti- 



mated velocity-dispersion errors are in excess of 15% of 
their maximum-likelihood velocity-dispersion values, we 
obtain 7 Z = —0.58 ± 0.15, which again differs only in- 
significantly from the value we find when including all 
lenses. Table Q] reports all parameters from this fit. 

4.2. Photometric parameter dependence 

From Figure [3l we see no significant overall redshift 
variation in the combined sample in either mean stel- 
lar mass M+, mean effective radius i? e , or mean Sersic 
index n. The lack of an apparent stellar-mass trend 
can be considered i n terms of t he lu minosity- function 
evolution results of iWake et al.l ([2006D . which show no 
significant evolution in the number density of massive 
elliptical galaxies from redshift z ~ 0.6 to the present 
epoch, once passive stellar evolution is taken into ac- 
count. This suggests that we are probing a similar range 
of galaxies across the redshift range spanned by our lens 
sample. The lack of a detected trend in effective ra- 
dius is also consistent with the modest size evolution 
seen for massive elliptica l galaxies at lower redshifts (e.g., 
Ivan der Wei et al1l2008| ) , again suggesting that our sam- 
ple represents a broadly similar population of galaxies 
across redshift. Of course, the question of merging com- 
plicates the identification of corresponding galaxy popu- 
lations across redshift, and we will return to this issue in 

m 
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Fig. 3. — Redshift variation of other observable quantities in the sample of 79 early-type SLACS (black diamonds) and BELLS (blue 
squares) lens galaxies. From upper left to lower right, these are: stellar velocity dispersion within one effective radius, estimated stellar 
mass, effective radius, Sersic index, ratio of lensing Einstein radius to lens-galaxy effective radius, and Einstein radius in kiloparsecs at the 
redshift of the lens galaxy. The solid black line shows the best-fit relation, and gray lines indicate the "1-sigma" error in the slope and 
zero-point of this relation. For all plots except the velocity-dispersion plot, the fit errors have been estimated via bootstrap resampling. 
For the velocity-dispersion plot, values have been scaled from the observational aperture (SDSS or BOSS) to one effective radius using the 
empirical relation of [Cappcllari ct al. (200(f), error bars indicate A\ 2 = 1, and dashed lines indicate the best-fit intrinsic RMS population 
scatter. Dot-dashed lines bound the region 2.2 < log 10 <r e < 2.5, and lenses with maximum-likelihood values outside this range are marked 
with red crosses. The velocity-dispersion evolution fit is done using the full X 2 (°"c) function for each galaxy, as described in the text. 
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We may also address the issue of heterogeneous sizes, 
stellar masses, and Sersic indices in our sample head- 
on, by introducing them as second independent vari- 
ables along with redshift in our fit for the dependence 
of the mass-profile parameter 7, since (in contrast to 
the errors on velocity dispersion) the statistical errors 
on M^, R e , and n are very small, and much more un- 
correlated with the dominant uncertainty in 7. For 
these second-independent-variable analyses, we use HST 
/-band de Vaucouleurs model m agnitudes and effective 
radii from iBolton et all ()2008| ) and Paper I. bersic 
indices are obtained by fitting two-dimensional PSF- 
convolvcd elliptical Sersic models to the same image 
data as used for the de Vaucouleurs photometry in those 
works, covering a 51" x 51" square region centered on the 
lens galaxies. Our stellar mass estimates are obtained by 
sc aling the passiv e ly-evo lving stellar-population model 
of iMaraston et aD (|2009f ) to the observed and Galactic 
extinction-corrected /-band magnitudes. 

We consider bi-variate population models of the form 

p(j\z,x; 70, 7z:7s,s 7 ) = J- x exp{-[7 - 

V Z7TS 7 

(70 + lz{z - 0.25) + lx {x - c x ))] 2 /{2s 2 )}, (8) 

where x = log lo (M*/M ), log 10 (/? G /l kpc), or log 10 (n) 
depending upon the second independent variable under 
consideration, with the constant c x = 11.7, 0.7, and 
0.7 for the three cases respectively. We grid the four- 
dimensional parameter spaces for these three cases in or- 
der to map out likelihoods and posterior probabilities, 
and show the results for the joint dependence parame- 
ters in Figure |U In all three cases, the dependence of 
7 on redshift is more significant than the dependence on 
the second variable. We evaluate the statistical improve- 
ments represented by these various two- variable models 
over the redshift-only model via the likelihood ratio test, 
expressed by the following A\ 2 statistic: 

Ax^ = 21n{max[£(7o,7*,7 a:) s 7 |{d i })]} 

-21n{max[£(7o,7„s 7 |{d 4 })]}. (9) 

The greatest improvement comes from the introduction 
of stellar mass as a second independent variable, for 
which Ax x = 5.07. The introduction of effective radius 
as a second dependent variable gives a fairly comparable 
Ax x = 4.65, while the introduction of dependence upon 
Sersic index gives a relatively insignificant Ax x = 1.82. 
Motivated by the evidence for stellar-mass dependence, 
we also compare a model with dependence on stellar mass 
alone to the initial model with dependence on redshift 
alone, and find the rcdshift-dcpcndcnce model preferred 
by Ax 2 — 7.16. We note also that our stellar mass es- 
timates are subject to much greater systematic uncer- 
tainty than our redshift estimates. To summarize, we 
see some evidence for shallower profiles in larger and/or 
more massive galaxies, but this cannot explain the evo- 
lutionary signal that we detect. 

4.3. Lensing aperture effects 

Figure [3] shows a non- negligible trend of increase with 
redshift in the Einstein radius, both in absolute units and 
relative to the effective radii of the lens galaxies. This 
trend is chiefly a consequence of the cosmic geometry 



of strong lensing for increasing lens redshifts: to leading 
order, a given galaxy has a characteristic angular scale 
for strong lensing that does not diminish with increasing 
redshift, whereas the apparent angular scale of a given 
physical length unit does diminish with increasing red- 
shift. Although the Einstein radius in kiloparsecs does 
not represent a physical property of the lens galaxy, we 
must consider the possibility that this variation of physi- 
cal Einstein radius with redshift impacts our mass-profile 
evolution measurement. If we assume as an alternative 
hypothesis that lens galaxies do not evolve structurally, 
our apparent evolution result must imply that the av- 
erage mass-density profile is shallower with increasing 
radii out from the center of the galaxy. At face value, this 
would be at variance with both a de Vaucou leurs or Sersic 
stellar profile and a iNavarro et al.1 (|1996| ) dark-matter 
profile, since both profiles become steeper at increasing 
radii. An intriguing possibility is that we might be see- 
ing direct evidence for an inflection zone between the 
stellar-dominated inner regions and the halo-dominated 
outer regions of our lens galaxies, which in certain com- 
bined baryon+halo models can exhibit a localized trend 
towards a shall ower profile with in creasing radius (see, 
e.g., Figure 8 of lGavazzi et al.ll2007t ). However, it should 
be noted that our lensing+dynamics analysis constrains 
an average density profile interior to the Einstein ra- 
dius, rather than a boundary value at the Einstein ra- 
dius, so this interpretation is not straightforward. The 
most direct means of addressing this question further 
would be through spatially resolved stellar kinematics 
of the SLACS and BELLS lenses, which in combination 
with the strong-lensing mass normalizations could test 
the redshift and aperture dependences separately. 

Contributions from large-scale structure fluctuations 
along the line of sight (i.e., beyond the host dark- matter 
halo scale of the lens galaxies) will not contribute signif- 
icantly to our results. The influence of large-scale struc- 
ture variance will increase with the redshift of the lensed 
background galaxy ZbGj and hence the contribution of 
this effect in SLACS and BELLS lenses with source red- 
shifts zbg is 1-0 will be significantly smaller than for 
lensed quasars or cluster-lensed giant arcs at higher red- 
shifts. Large-scale structure fluctuations in the lensing 
convergence k (projected surfa ce density scaled to the 
critic al density for lensing: e.g.. iNaravan fc Bartelmannl 
119961 ) are expect ed to be below the 1% level over this red- 
shift range (e.g. JDalal et al.ir2005[ ). Hence these fluctua- 
tions are sub-dominant to the lensing Einstein measure- 
ment errors, which are in turn subdominant to velocity- 
dispersion uncertainties. In addition, since these will be 
fluctuations about the mean density, they are a source 
of noise rather than bias in our measurement, and will 
average out with increasing sample size. 

4.4. Spectroscopic and lensing selection effects 

lArneson et al.l (|2012[ ) have quantified the combined 
bias effects of spectroscopic candidate selection and 
strong-lens confirmation as a function of mass-density 
profile parameters for lens surveys with the charac- 
teristics of SLACS and BELLS, using detailed Monte 
Carlo simulations. The SLACS and BELLS simula- 
tions are distinguished primarily by differences in spec- 
troscopic emission-line detection depth (BELLS com- 
ing from deeper BOSS data), background-galaxy red- 




A 

V 
_<D 

£ 

D 

</) 



D 



"O 



Q. 



0.03 
0.02 
0.01 
0.00 
-0.01 

-0.02 
-0.03 



1 1 1 


. , ........ . | 


r o <> 


□ 

□ ; 






o«>o % 


u □ 5 — : 
o % D □ ; 


□ 

' 1 ' 


a i 



Fig. 4. — Posterior probability contours enclosing credible regions for the joint dependence of the population mean logarithmic mass 
profile slope 7 on: redshift and stellar mass (left) ; redshift and e f fectiv e radius (center); and redshift and Sersic index (right). Stellar 
masses are estimated by scaling passively evolving Maraston et al. (2009) models to the observed HST F814W magnitudes. Dashed lines 
indicate the lines of no dependence of the population mean 7 value on redshift (horizontal) or on the various second independent variables 
(vertical). Zero-point and intrinsic population scatter parameters have been marginalized in all cases. 

shift (BELLS being higher on average), and spectroscopic 
fiber size (BELLS having the 2"-diameter BOSS fiber, 
SLACS having the 3"-diamctcr SDSS fiber). Figure 5 of 
that work presents the relative selection and confirmation 
probabilities for lenses as a function of angular Einstein 
radius 6*e and mass-profile slope 7 for both SLACS-likc 
and BELLS-like surveys. The 6e dependence incorpo- 
rates the dependence on velocity dispersion, lens redshift, 
and source redshift in the particular scaled angular com- 
bination that is directly relevant for strong lensing. We 
use the data from this work to determine the magnitude 
of any possible lensing selection bias in our mass-profile 
evolution measurement. 

For each len s in our sample , we se lect all simulated sys- 
tems from the lArneson et ail (|2012f) data that are within 
±0.1 dex in 6*e of that lens's 6*e value, and that meet 
the simulated criteria of being both spectroscopically se- 
lected and modelable as strong lenses. Since the parent 
population of simulated lenses is distributed uniformly in 
7 from 1 to 3 (the full range of mathematically allowable 
values), the histogram of these selected subsets defines 
the relative probability of a simulated lens appearing in 
our observational sample as a function of 7. We then 
assume an intrinsic distribution of 7 values given by the 
best-fit population parameters from $3] above, evaluated 
for the redshift of the real lens under consideration. We 
weight this intrinsic distribution by the 7-dependent se- 
lection probability, and compute the mean 7 value over 
the selected subset of simulated lenses with this combined 
weighting. Finally, we subtract off the intrinsic mean 7 
value to determine the expected bias in 7 for the real 
lens under consideration. Figure [5] shows the results of 
this calculation for all lenses in our observational sample. 
The expected selection bias in 7 is at the level of ~ 0.01, 
well below the evolution signal that we detect, and shows 
no significant redshift dependence or variation between 
the SLACS and BELLS samples. 



5. DISCUSSION AND CONCLUSIONS 

Comparing our m ass-profile evolution result to that 
of lRuff et all (l20Tl . we find that our best value for j z 
is more than twice as large. However, the significance 
of the discrepancy, given the combined uncertainties in 
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Fig. 5. — Expected mass-profile parameter biases for SLACS 
(black diamond) and BELLS (blue square) lenses, based on com- 
bined spe ctroscopic selection a nd lens confirmation biases as sim- 
ulated bv lArneson et al.l J2012I ). 

both measurements, is less than 2-sigma. Furthermore, 
we note that the magnitude of the mass-profile evolution 
detected in this work, although expressed per unit red- 
shift , is detected over a redshift interval less than one PI 
The lRuff et al.l (|2011l) analysis includes lens systems out 
to redshift z « 1 from the Lenses Structure and Dy- 
namics Survey (LSD: e.g . , iTreu fc Koopmansl[200l [2004 : 
iKoopmans fc Tret] 120031 ). Hence if the bulk of evolution 
is concentrated towards the more recent half of the inter- 
val since redshift unity (which encompasses 65% of the 
cosmic time elapsed), the two results would be brought 
into further agreement. We also note that our result is 
consistent with t h e lack of significant evolution found in 
IKoopmans et al.l (|2006| ) given the uncertainties in that 
work. Finally, we note that the tension between the two 
measurements is somewhat lessened by the possibility of 
an additional mode of variation of 7 with lens-galaxy 
mass and/or size. 

13 Az 0.6 measured from the minimum SLACS lens redshift 
of 0.06 to the maximum BELLS lens redshift of 0.66, or Az ~ 0.3 
measured from the median SLACS redshift of 0.2 to the median 
BELLS redshift of 0.5. 
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One possible explanation for evolution towards steeper 
mass profiles with cosmic time could be an ongoing 
contribution from dissipative processes, as suggested by 
iRuff et aTl (|2011f ). However, baryonic dissipation suf- 
ficient to significantly modify the global mass-density 
structure would necessarily be accompanied by appre- 
ciable star-formation, which is at variance with the uni- 
formly old stellar populations that are characteristic of 
massive early-type g alaxies (e.g., I Thomas et al.l 120051 : 
iMaraston et al.lT2009l ). An alternative possibility is that 
hierarchical dry-merging processes can lead to an evo- 
lution in the inner mass-density profile. Cosmological 
simulations of the formation and assembly of galaxies 
suggest that massive (M* > 10 11 M Q ) ellipticals will typ- 
ically have assembled only h alf of their z = s tellar 
mass at a redshift of z = 0.8 (|De Lucia et al.l [20061 ). On 
smaller-scales that resolve th e density stru cture within 
galaxies, bot h analytic work (IDehnenl 120051) and numer- 
ical studies (jFulton fe Barnes! 1200 It iKazantzidis et al.l 
I2006D find that halo core profiles are preserv ed under ma- 
j or me rging. However, in the simulations of iNipoti et al.l 
(|2009D . which focused on strong- lensing observables in 
the evolution of the central mass-density profile un- 
der the merging of isothermal progenitors, the specific 
channel of off-axis major mergers produces a system- 
atic evolution towards steeper-than-isothcrmal remnant 
profiles, with a magnitude comparable to the integrated 
evolution seen here in our SLACS+BELLS lens sam- 
ple. A possible scenario that emerges is hence one in 
which the evolution in the distribution of stellar trac- 
ers (and hence in the observed mass-size relation) in 
massive early- typ e galaxies is driven primarily by minor 
dry mergers (e.g.. iNaab et al.|[200llHopkins et alj|2009t 
Ivan Dokkum et al .1 12010D at higher redshifts, while an 
evolution in their total mass-density structure is driven 
by major dry mergers at lower redshifts, with off-axis 
merger configurations playing the most important role. 
An alternative explanation could be an evolution to a 
higher fraction of sat ellite (relative to central) galaxies at 
lower redshifts (e.g., IWake et all 120081 : I Wetzel fc White! 
l2010f ). with consequently increased central mass concen- 
trations from the effects of tidal stripping. 

When galaxy mergers are considered, the question of 
matching present-day galaxies with their population pro- 
genitors at higher redshift becomes more complicated. 
Our result should therefore be regarded as a constraint 
on the recent structural evolution of galaxies within the 
approximate range 10 11 < M±/Mq < 10 12 , with the un- 
derstanding that if major galaxy mergers are a significant 
driver of that evolution, the identity of a given galaxy is 
not fixed by its stellar mass even in the absence of star 
formation. A single major dry merger with a 1:1 mass 
ratio will place the product galaxy 0.3 dex higher in stel- 
lar mass than its two progenitor galaxies, which repre- 
sents a relatively modest change compared to the range 
of masses spanned by galaxies in the universe. If we fur- 
thermore assume that similar physical processes are at 
work for galaxies just beyond the reach of our sample 
in stellar mass, we can be confident that our analysis is 
generally applicable to the structural evolution of mas- 
sive galaxies. With more detailed observations, a more 
controlled sample, and a tighter quantitative connection 
to galaxy merger simulations, we can make our current 
result more precise as to the implied rate and nature of 



mergers. 

We have many avenues to expand our observational 
study of the density-structure evolution of massive el- 
liptical galaxies using the combined SLACS and BELLS 
samples. First, we will pursue an analysis based upon 
more physically motiv ated star-plus-halo model s of the 
lens galaxies (see, e.g.. Uiang fe Kochanek|[2007l ). rather 
than the somewhat ad hoc class of power-law models 
considered here. Second, we will pursue deeper ground- 
based spectroscopy of the BELLS lens sample, since the 
BELLS velocity-dispersion uncertainties are the domi- 
nant limitation to the statistical precision of our cur- 
rent result. Third, deep and high-resolution multi-band 
imaging of BELLS lenses can be obtained, since cur- 
rently our only multi-band imaging data for the BELLS 
lenses is from the SDSS, and is affected by large pho- 
tometric errors and background-galaxy flux confusion. 
These deeper data can be used to model the lens-galaxy 
stellar populations in greater detail, and to determine 
stellar masses that are more accurate than the simple 
luminosity-scaling estimates that we have adopted in this 
paper. Better stellar-mass estimates would also allow for 
a meaningful separation of the stellar and dark-matter 
contributions to the total lensing mass. Fourth, we can 
incorporate mass-profile constraining information from 
the resolved lensed images in addition to stellar dynam- 
ics. Finally, the availability of new BELLS lens can- 
didates from continuing BOSS observations should lead 
to a ten-fold increase in the sample of lenses at higher 
redshift, allowing us to more robustly measure the joint 
dependence of early-type galaxy structure on the two in- 
dependent variables of mass and redshift. 
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